Micromotion-induced Limit to Atom-Ion Sympathetic Cooling in Paul Traps 
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We present and derive analytic expressions for a fundamental limit to the sympathetic cooling 
of ions in radio-frequency traps using cold atoms. The limit arises from the work done by the 
trap electric field during a long-range ion-atom collision and applies even to cooling by a zero- 
temperature atomic gas in a perfectly compensated trap. We conclude that in current experimental 
implementations this collisional heating prevents access to the regimes of single-partial-wave atom- 
ion interaction or quantized ion motion. We determine conditions on the atom-ion mass ratio and 
on the trap parameters for reaching the s-wave collision regime and the trap ground state. 
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The combination of cold trapped ions and atoms [IH^I 
constitutes an emerging field that offers hitherto unex- 
plored possibilities for the study of quantum gases. New 
proposed phenomena and tools include sympathetic cool- 
ing to ultracold temperatures 0, , charge transport in 
a cold atomic gas [ll], [HI], dressed ion-atom states [l3r- 
llH ]. local high-resolution probes [13, EH and ion-atom 
quantum gates (lil. |20|. 

In contrast to atom traps based on conservative forces, 
Paul traps employ radiofrequency (RF) electric fields to 
create a time-averaged secular trapping potential for the 
ion [2l| . The time- varying field can pump energy into the 
system if the ion's driven motion is disturbed, e.g., by a 
collision with an atom [22j. The kinetic behavior of an 
ion in a neutral buffer gas has been observed in numerous 
experiments [23H3Q| ■ The ion's equilibrium energy distri- 
bution was predicted analyticall y 1311. as well as using 
Monte-Carlo techniques [111 [27, 32l 33|. and recently a 
quantum mechanical analysis has been performed [34j . 
For atom-ion mass ratios below a critical value, the ion 
is predicted to acquire a stationary non-thermal energy 
distribution with a characteristic width set by the coolant 
temperature [33j]. 

The RF field drives micromotion of the ion at the RF 
frequency. At any position and time, the ion's velocity 
can be decomposed into the micromotion velocity and 
the remaining velocity of the secular motion. Consider 
the simple case of a sudden collision with an atom that 
brings the ion to rest, a process which in a conservative 
trap would remove all kinetic energy. Immediately af- 
ter such a collision the ion's secular velocity is equal and 
opposite to the micromotion velocity at the time of the 
collision. This implies that the ion's secular motion can 
be increased even in a collision that brings it momentar- 
ily to rest. Hence for cooling by an ultracold atomic gas 
0,01, the energy scale of the problem is no longer set by 
the atoms' temperature but by the residual RF motion 
of the ion, caused, e.g., by phase errors of the rf drive 
or by dc electric fields which displace the trap minimum 
from the rf node [HI]. Such technical imperfections have 



limited sympathetic cooling so far, with the lowest in- 
ferred ion temperature on the order of 0.5 mK 0. Full 
quantum control in these systems 19|, 13; 011 the other 
hand, likely requires access to the smaller temperature 
scales ftw/kB ~ 50^iK for the trap zero-point motion and 



E s /kB ~ 50 nK for the s-wave collision threshold 

In this Letter, we show that even with the atomic 
gas at zero temperature and in a perfectly dc- and rf- 
compcnsated Paul trap, a fundamental limit to sympa- 
thetic atom-ion cooling arises from the electric field of the 
atom when polarized by the ion, or equivalently, the long- 
range ion-atom interaction. As the atom approaches, it 
displaces the ion from the RF node leading to micromo- 
tion, whose interruption causes heating. A second non- 
conservative process arises from the non-adiabatic mo- 
tion of the ion relative to the RF field due to the long- 
range atom-ion potential; during this time, the trap can 
do work on the ion and increase its total energy. We 
find that, in realistic traps, the work done by the RF 
field dominates the effect of the sudden interruption of 
the ion's micromotion and leads to an equilibrium energy 
scale that, for all but the lightest atoms and heaviest 
ions, substantially exceeds both the s-wave threshold E s 
and the trap vibration energy hw. Our analysis shows 
that current atom-ion experiments [H-0| will be confined 
to the regimes of multiple partial waves and vibration 
quanta, and indicates how to choose particle masses and 
trap parameters in order to achieve full quantum control 
in future experiments. Our analytical results are sup- 
ported by numerical calculations that furthermore reveal 
that in those collisions where the RF field removes energy 
from the system, the atom becomes loosely bound to the 
ion, leading to multiple subsequent close-range collisions 
until enough energy is absorbed from the RF field to eject 
the atom and heat the ion. 

We consider a classical model and later confirm that 
this assumption is self-consistent, i.e. that the energies 
obtained from the model are consistent with a classical 
description. An atom of mass m a approaches from infin- 
ity to an ion of mass m, held stationary in the center of 
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an RF quadrupole trap. At sufficiently low collision en- 
ergies, the angular-momentum barrier will be located far 
from the collision point, and, once it is passed, the colli- 
sion trajectory will be nearly a straight line. We initially 
assume this trajectory to be along an eigenaxis of the RF 
trap, resulting in a true ID collision in a total potential 
given by V (ri,r a ,t) — e£ (fj, t) ri/2+U(r), where e is the 
ion's charge, rj and r a are the ion and atom locations, 
respectively, r = 7*j — r a is the ion-atom distance, and 
£ (r.i , t) — gri cos (Sit + <f>) is the RF electric field of the 
ion trap at frequency SI, parameterized by its quadrupole 
strength g. The ion- atom interaction potential at large 
distances is given by U(r) — — C4/(2r 4 ) [ll| with C4 the 
atom's polarizability, and modeled as a hard-core repul- 
sion at some small distance. Since in a three-dimensional 
collision the atom can approach the ion along directions 
that are perpendicular to the RF field or where the RF 
field vanishes, we expect our analysis to overestimate the 
heating by a factor of order unity, which we address be- 
low. 

As the ion is pulled from the trap origin by the long- 
range interaction U with the approaching atom, the os- 
cillating electric field £ causes it to execute sinusoidal mi- 
cromotion with amplitude (jTj/2, where q = 2eg / (rriiSl 2 ) 
is the unitless trap Mathieu parameter. As long as the 
motion of the ion relative to the atom remains slow rel- 
ative to the average RF micromotion velocity o Mm , the 
ion's equations of motion during one RF cycle will remain 
linear, and the secular motion of the ion during each RF 
cycle can be described in terms of an effective conser- 
vative secular potential V s — \m i uj 2 r 2 + U(r), where 
u) rj qSl/2 3 / 2 is the secular frequency of the ion within 
the trapping potential. Associated with V s are the char- 
acteristic length scale R = C^ 6 / (m^w 2 ) 1 ^ 6 at which the 
interaction potential U is equal in magnitude to the trap 
harmonic potential, time scale T = 2ix /w, and energy 
scale E R = \ mi u 2 R 2 = ^m^^d) 1 / 3 . 

In collisions with light atoms, m a < m*, we expect 
the ion to stay confined close to the trap origin. In this 
case, the ion-atom distance r(t) is governed solely by the 
motion of the atom in the ion-atom potential U, 
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Figure 1: Trajectories of an ion r% (t) and an atom r a (t) dur- 
ing a classical one-dimensional low-energy collision. The atom 
of mass m a approaches the ion of mass rrn = 2m a held in the 
center of a RF trap with secular frequency cj — 2ir/T and 
Mathieu parameter q — 0.1, leading to a hard-core collision 
at n — r a — r c , t — and RF phase <j>. For <f> = n/2 (dotted 
lines), the trap field adds energy to the system, causing heat- 
ing. For (j> — 3tt/2 (solid lines), the RF field removes energy, 
binding the atom to the ion and causing further collisions at 
various RF phases until enough energy is accumulated to eject 
the atom. 
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Figure 2: Change AE in the total secular energy of a colliding 
ion-atom system (jn a jm,i — 1/2, q — 0.1), as a function of 
the RF phase cj) during the first hard-core collision. For < 
<f> < n, the system undergoes only one collision with energy 
gain comparable to the analytic prediction (Q (solid line). 
For 7r < (f> < 2tt, the atom is bound and undergoes several 
collisions with the ion before eventually escaping. 



where the collision occurs at time t = and fi — 
mim a /(mi + m a ) is the reduced mass of the system. The 
displacement r c of the ion from the center of the trap at 
the time of the hard-core collision with the atom can then 
be estimated by integrating the effect of the force 2Ci/r b 
exerted by the atom on the ion trapped in its secular po- 
tential, yielding r c rj 1.11 (m a /mi) 5 ^ 6 R. In collisions 
with heavy atoms (rrii < m a ) on the other hand, the 
ion responds quickly to minimize the total secular po- 
tential energy until, at (rj,r a ) = (0.29, 1.76)7?, the de- 
formed ion's equilibrium position becomes unstable and 
the light ion quickly falls towards the atom with the col- 



lision occurring at the ion displacement r c sa 1.76-R. In 
general, we may express the collision point as r c /R = 
r c (q, m a /m t ) / (1 + mi/m a ) 5/6 where 1 < r c < 2. 

Around the collision point r c , we may expect the en- 
ergy of the system to change as long as the ion moves 
non-adiabatically relative to the RF field, including the 
interval —t\ < < t\ around the collision at t — 0, dur- 
ing which the ion's velocity r, is greater than its average 
micromotion velocity sa ujr c . In this regime, the trap 
RF field can be thought of as a time-dependent pertur- 
bation to the ion-atom potential, doing work on the ion 
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equal to 



W 



8 (n (t) , t) ■ h (t) dt. 



(2) 



The change in the system's energy will depend on RF 
phase <p at the time of the close-range ion-atom collision, 
taking on the maximal value W max for <j> = <j) m ax- For 
fi 3> Vfim, we ma y neglect the effect of the electric field 
on the ion's trajectory and approximate the ion's position 
in terms of the free collision trajectory r(t), Eq. [TJ as 
Ti k r c — m a r /(mi + m a ). In this case, the work done by 
the RF field can be written as W as K sin 6. where 
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is the characteristic scale for the work done on the ion 
by the RF field. The maximal energy gain W max a K 
occurs for <fi = <fi max s» 7r/2, corresponding to the RF 
field changing sign at the time of the collision. The 
non-adiabatic condition |rj| > \v^ m \ is equivalent to 

3<?f2t < (2/f c ) 3 ^ 2 , which, for the practically relevant val- 
ues of q < 0.5, will always include the region \Qt\ < 0.8 
where the dominant contribution to the integral in ([3]) 
occurs. Consequently, we may extend the limits of in- 
tegration to ±oo to obtain K/Wo ~ 1.82f c — l.&'&q 1 ^, 
implying 0.7 < K/Wo < 2.8 for q < 0.5 and at all mass 
ratios. Under the same conditions, K is at least three 
times larger than the ion's average micromotion energy 



at the collision point 



rriiV,, 



l /2 and the gradual 



energy change (JSJ) dominates the effect of the sudden in- 
terruption of the ion's micromotion. Intuitively, at low 
collision energies, the ion-atom potential dominates for 
a longer time during which the RF field does more work 
on the system. Since the trap electric field must be in- 
creased for higher RF frequencies to preserve the ion sec- 
ular potential, the heating increases with a decrease in 
the Mathieu parameter (K > 15 x E^ m for q < 0.1). 

Since W corresponds to the difference in the work done 
by the RF field during the incoming and outgoing parts of 
the collision, the energy change will depend on the phase 
<\> of the RF field at the time of the hard-core collision. 
For < 4> < 7r, the RF field accelerates the collision 
partners towards each other, increasing the system en- 
ergy; for 7r < <fi < 2tt, the RF field opposes the collision, 
doing negative work and causing the atom to be bound 
to the ion with binding energy on the order of —Wo (Fig- 
ure [TJ). Since the r -4 potential does not possess stable 
orbits, bound ion-atom trajectories will include further 
close-range collisions. Depending on the RF phase dur- 
ing each subsequent collision, the system will gain or lose 
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Figure 3: Maximal energy gain Wmax in a low-energy ID ion- 
atom collision as a function of the atom/ion mass ratio m a /mi 
with Mathieu parameter q = 0.1 (a) and as a function of q 
with nia/riii — 1/2 (b). Data points are numeric calculations. 
The solid lines correspond to W ma x = K, while the dashed 
lines represent W ma x = Wo- 
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Figure 4: (a) Distribution of numerically computed secu- 
lar energy gains in 2883 random low-energy collisions be- 
tween a free 87 Rb atom and a 174 Yb + ion held in the three- 
dimensional trap from [2j|, in units of the 3D micromotion 
heating energy scale Wo D = 4 Wo / (3tt). A sample ion-atom 
collision trajectory is shown in the inset (b). 



energy on the order of Wo, leading to a random walk 
in energy until the atom finally unbinds. Then AE de- 
pends sensitively on the RF phases at each hard-core 
collision spaced in time by many RF cycles, leading to 
a sharp dependence of AE on the RF phase (j> during 
the initial collision for ir < <j> < 2ir (Figure^. Using 
the approach below, we calculate the net average energy 
gain for < 4> < 2n and m a /mi = 1/2, q = 0.1 below as 
(AE) = 1.0W . 

To quantitatively verify the above heating model, we 
numerically calculated classical trajectories of low-energy 
ID collisions as a function of (j>, m a /mi and q (Figures 
[5] and [3]) . The ion was initially held at the trap center 
while the atom followed the analytic trajectory ([TJ) . At 
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a critical ion-atom distance ro, the ion was displaced so 
that the trap's secular force would balance the atom's 
attraction, while keeping its velocity zero. Choosing 
r = 3.8 (1 + mi/rria) 1 ^ 3 R ensured that for all the trap 
parameters in this paper, the ion's intial micromotion 
energy was smaller than 10~ 5 Wo and both the atom's ki- 
netic energy and U (ro) were smaller than 1Q~ 2 Wq. The 
equations of motion were integrated using the Dormand- 
Prince explicit Runge-Kutta method: away from collision 
points, the motion was integrated as a function of time 
while near the collisions, the ion-atom distance r was 
used with a hard-core radius e = 10~ 3 (1 + mi/rria) 1 ^ R. 
The integration was stopped when the atom reached a 
distance r a = 2.1ro much larger than the ion motion, at 
which point the total secular energy of the system was 
evaluated. We confirmed the accuracy of our integration 
by replacing the RF potential with a time-independent 
secular potential and confirming energy conservation at 
the level of 10~ 3 W . 

Figures [3^ and [3t> show the numerically calculated 
maximal energy gain W max in the initial collision as 
a function of m a /m,i and q. For q < 0.1, the cal- 
culated W ma x is within 30% of the analytic prediction 
W max = K: for q = 0.5, m a /mi = 1/2, W max increases 
to 1.7 K, partly because the analytic result does not in- 
clude the micromotion interruption. The heating is in- 
sensitive to the the hard-core radius: a tenfold increase or 
decrease in e changes W max by less than 1%. Our results 
are also very robust with respect to the initial conditions: 
the average energy gain in Figure [2] changes by less than 
5% if the ion is started at rest in the center of the trap 
and the atom at rest a distance R away. 

In three-dimensional linear quadrupole RF traps, the 
electric field £ and the ion position rj in © are vectorial 
quantities. Assuming that the collision trajectory is still 
nearly one-dimensional, averaging over the atom's ap- 
proach direction rescales the work done by the RF field, 
Eq. 12]), by 4/ (3ir), leading to a natural 3D heating scale 
Wq D = 4 Wo/ (37r). To check this, numerical simulations 
were done for the three-dimensional quadrupole RF trap 
from Ref. 0. We considered cold 87 Rb atoms that have 
passed the angular momentum barrier at a large distance 
and are colliding head-on with a 174 Yb + ion from a ran- 
dom direction and at a random time. The initial condi- 
tions for numeric integration were chosen by analogy to 
the ID case. Figure |4Jd shows a sample ion-atom trajec- 
tory including multiple collisions. Due to a difference in 
the axial and radial frequencies of the ion trap, the colli- 
sion trajectory precesses about the y axis, while remain- 
ing nearly one-dimensional close to the collision points. 
A histogram of the final system energies E after the atom 
is ejected to infinity is shown on Figure 0^, with an aver- 
age energy gain of 0.9Wq D , confirming that the heating 
is similar in one and three dimensions. 

Table Q] shows Wq together with the s-wave threshold 
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174 Yb+ 87 Rb [2] 


200 


0.013 


44 


4.6 


9.6 


540 


87 Rb+ 87 Rb [7] 


350 


0.24 


79 


7.7 


17 


210 


138 Ba+ 87 Rb [4] 


200 


0.11 


52 


4.5 


9.6 


150 


40 Ca+ 87 Rb [6] 


110 


0.20 


200 


2.6 


5.3 


50 


i74 Yb+ i72 Yb fl] 


67 


0.14 


44 


1.6 


3.2 


41 


174 Yb + 40 Ca [gj 


250 


0.25 


270 


14 


12 


32 


174 Yb + 23 Na 


50 


0.30 


710 


4.7 


2.4 


1.5 


174 Yb+ 7 Li 


50 


0.30 


6400 


14 


2.4 


0.24 



Table I: The quantum s-wave energy limit {E s ), the corre- 
sponding radial DC electric field E S DC at which the ion micro- 
motion energy is equal to E s , the ion trap vibrational quan- 
tum (ftw), and the micromotion- induced energy scale Wo in 
various ion-atom systems. The Yb + +Rb and Ba + +Rb sys- 
tems have a sufficiently large ratio of elastic to inelastic colli- 
sions to permit cooling [3, i|. In the Yb + +Yb and Rb++Rb 
systems, it may be possible to make charge exchange colli- 
sions endothermic by choosing the proper isotope combina- 
tion. The Yb + +Ca system exhibits a very large inelastic col- 
lision rate. 



energy E s = h 4 / (2p?Ca) and the energy Huj of a trap vi- 
brational quantum in various experimental systems. In 
the current systems P, 0, 0, H, 0| , the atom number per 
characteristic volume R 3 is much smaller than one, while 
the atom kinetic energy is much smaller than Wo, and 
we expect binary ion-atom collisions to be well described 
by our model. Given the typical duration of a collision of 
about T/2 and typical ion-atom collision rate coefficients 
0, corrections due to three-body collisions are expected 



at critical atom densities on the order of n c 



10 



14 



At these densities, three-body inelastic processes become 
significant, leading to strong ion heating and loss as was 
recently observed 0. Thus best sympathetic cooling is 
expected in the density regime n < n c , where it is limited 
by the two-body energy scale Wo derived above (Eq. [3} 
to temperatures T m i n ~ Wo/ks- In currently realized 
systems, Wo is more than one order of magnitude larger 
than the trap vibrational quantum and almost three or- 
ders of magnitude larger than the s-wave scattering limit. 

Since, for light atoms, the ratio of of heat- 
ing to the s-wave collision threshold scales as 
W)/fiw oc (wCi) 4 ^ 3 rria 1 /irii and the ratio of heating 
to the trap vibration quantum scales as Wa/(hu) oc 
(ujC/l) 1 ^ 3 m a /m.i, our model predicts that micromotion 
heating could be mitigated using light atoms and heavy 
ions trapped in weak RF traps, limited by the control of 
DC fields (since E S DC oc w). In particular, with control 
over the DC electric fields on the order of 10 mV/m - an 
order of magnitude better than current state-of-the-art 
[7| - the Yb + /Li system may enter the s-wave regime. 
Heating could also be decreased by employing Raman 
transitions to produce molecular ions that are sufficiently 
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bound so as not to be affected by micromotion [34j. An- 
other option may be the use of a n op tical trap for the 
ion, as was recently demonstrated |3a |. 
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